Preparation of Decoherence Free Cluster States with Optical Superlattices 
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We present a protocol to prepare decoherence free cluster states using ultracold atoms loaded in 
a two dimensional superlattice. The superlattice geometry leads to an array of 2 x 2 plaquettes, 
each of them holding four spin- 1/2 particles that can be used for encoding a single logical qubit 
in the two-fold singlet subspace, insensitive to uniform magnetic field fluctuations in any direction. 
Dynamical manipulation of the supperlattice yields distinct inter and intra plaquette interactions 
and permits to realize one qubit and two qubit gates with high fidelity, leading to the generation 
of universal cluster states for measurement based quantum computation. Our proposal based on 
inter and intra plaquette interactions also opens the path to study polymerized Hamiltonians which 
support ground states describing arbitrary quantum circuits. 



I. INTRODUCTION 

Quantum technology, in particular quantum informa- 
tion processing and quantum metrology, requires the pre- 
cise preparation of quantum states that outperform a 
given task better than any classical strategy. As shown in 
recent years, the unprecedent control and precision pro- 
vided by ultracold gases in optical lattices makes these 
systems optimal candidates for such a technology. Fur- 
thermore, since the dynamical control over the optical 
lattice parameters permits the simultaneous coupling be- 
tween nearest atom-lattice sites, these systems are also 
increasingly used as quantum simulators to mimic dis- 
tinct complex condensed matter Hamiltonians [T]. 

The controlled generation of double well lattices, i.e. 
lattices whose unit cells contains two sites, has opened 
the possibility to isolate and address individually pairs 
of atoms, and hence to manipulate the interactions be- 
tween them. Seminal results are the demonstration of 
controlled exchange interaction between pairs of neutral 
atoms in an optical lattice when the atoms are forced 
to be in the same location [5], and the demonstration of 
superexchange interactions [5], showing that the inter- 
actions between atoms trapped in two adjacent sites of 
the optical lattice can be made analogous to the interac- 
tions between atomic spins in magnetic materials. While 
Ref. [2] sets a basis to perform in a controlled manner 
two qubit gates between neighboring atoms in the dou- 
ble well lattice, Ref. [3] opens a direct path towards the 
realization of low-temperature quantum magnets and a 
variety of many-body spin models with ultracold atoms. 

The atomic interaction control achieved with optical 
lattices has also direct applications to quantum compu- 
tation. A particularly well suited approach that exploits 
the innate massive parallelism of such systems to per- 
form quantum computation is the Measurement Based 
Quantum Computation (MBQC), where information is 



processed by means of a sequences of measurements on a 
highly entangled initial state. It requires the capability 
to create universal cluster state, that is, a multipartite 
quantum state able to reproduce any entangled quantum 
state in 2D, and to perform local single qubits opera- 
tions. Using as qubits two internal states of atoms in a 
2D optical lattice, it is possible to create a highly entan- 
gled quantum state by means of controlled collisions [3], 
which is indeed a prerequisite for the generation of uni- 
versal cluster state. 

Although neutral atoms couple weakly to the environ- 
ment and they have relatively long coherence times com- 
pared with the time scale associated with the achievable 
coupling strength, when atoms are brought to an entan- 
gled state decoherence will rapidly destroy any quantum 
superposition of atoms. The larger the entangled system 
is, the faster it will decohere. To fight against decoher- 
ence one should prepare the atoms in quantum states that 
are robust against external perturbations. For periodic 
arrays of double- wells [3 [6l |7] , resilient encoding schemes 
using two internal states (two Zeeman levels) of the atoms 
have been proposed. In these schemes, each double-well 
traps two two-level particles to encode a logical qubit. 
The logical space is spanned by the singlet and triplet 
states of the two spin 1/2 particle along the quantization 
axis (here denoted by z) = ^ (|t) U) - |i) |T)) and 

\To) = 73 (IT) U) + li) IT))}- Since these states have a 
zero z-component of the total spin, such encoding is in- 
sensitive to fluctuations of the magnetic field along the 
quantization axis. In practice, such an encoding scheme 
is very well suited for robust controlled interactions along 
one, lets say horizontal, direction. To create a universal 
2D cluster states (cluster states in ID are not universal), 
interactions and hence entanglement between neighbor- 
ing atoms along both, the horizontal and the vertical 
directions, should be performed in such a way that ro- 
bustness is preserved. In the above encoding, interac- 



2 



tions along the vertical direction, will leave the subspace 
of zero spin component along the quantization axis, be- 
coming very fragile in front of external magnetic field 
fluctuations noise. In this paper we show how this limi- 
tation can be overcome by using 2D optical superlattices. 

In passing, let us point out that superimposing sec- 
ondary optical lattices (or " superlattices" ) on top of the 
primary ones to further modify the potential in which 
the atoms are trapped permits in general to create poly- 
merized lattices. By polymerized lattices we mean lat- 
tices consisting of weakly coupled groups of neighboring 
atomic sites denoted as plaquettes. An example of a poly- 
merized lattice is a square lattice made of smaller squares. 
The intra-plaquette interactions in such lattices might be 
strong and may even be designed to include many(> 2)- 
body terms, while the inter-plaquette interactions might 
be much weaker. Polymerized lattices allow for instance 
to engineer Valence Bond Solids on demand, to study 
topological spin liquids and one might envisage them as 
potential quantum circuits. It is also very appealing to 
try to use the plaquettes as qubits or qudits (elementary 
systems with more than two internal states) for quantum 
information processing and implement quantum logical 
gates, quantum protocols, and quantum error correction 
in such systems by employing either interatomic inter- 
actions or/and interactions with external (electric, mag- 
netic, laser) fields. 

Here we take advantage of the 2 dimensional super- 
lattices to present new schemes to prepare universal 2D 
cluster states using the plaquettes as logical qubits. The 
superlattices creates a periodic array of plaquettes, i.e., 
2x2 potential wells (as shown in Fig. [T|), each of them 
filled with an atom with two internal degrees of freedom 
(spin 1/2 particle). On each plaquette, we encode a sin- 
gle logical qubit using the two-fold singlet subspace of the 
four 1/2-spins, as shown in Fig. [l] Thus, by doubling 
the physical resources in comparison with the 2 physical 
qubit encoding previously mentioned, we obtain the de- 
sired encoding that is decoherence free against uniform 
magnetic field fluctuations in arbitrary directions. 

The encoding scheme using the singlet subspace of four 
qubits has been previously studied for the quantum dot 
systems, and it is also called the "supercoherent qubit" 
[SI m [To] . Also, it has been shown that in such a config- 
uration, tunable Heisenberg superexchange interactions 
(between neighboring spins, including the diagonal and 
off diagonal ones) are sufficient for universal quantum 
computation [Il][12j. In the model we propose here, the 
generation of a universal cluster state demands: (i) the 
ability to perform one qubit gates to prepare all logi- 
cal qubits in the initial state \+) = ^/^/{2)[\0) + |1)], 
and (ii) the realization of controUed-phase gates, U = 
diag(l, 1, 1, —1), between nearest logical qubits, i.e. be- 
tween plaquettes, to create a maximally entangled 2D 
cluster state. Notice that the superlattice geometry of 
Fig. [l] does not induce interactions along the diagonals 
sites on each plaquette, and thus, the most challenging 
ingredient for universality is, indeed, the two qubit (two 



plaquette) gate. In this article, we propose three dif- 
ferent approaches to couple the logical qubits that ei- 
ther preserve the singlet subspace at the end of the gate 
operation or keep the state within the singlet subspace 
even during the whole completion of the gate. In all ap- 
proaches, we take into account realistic available tools 
and discuss the practical limitations in optical superlat- 
tices. In our flrst approach, we exploit the additional 
vibrational mode of the optical trap to facilitate the logi- 
cal coupling gate. In our second approach, we extend the 
earlier proposals [HI [12] by removing the requirement of 
equal coupling strengths for all six pairs within the pla- 
quette (more feasible for 2D optical superlattices), while 
we still obtain the effective Hamiltonian within the log- 
ical subspace sufficient for universal gates. In our last 
approach, we include tunable Ising-type interactions be- 
tween neighboring spins (attainable with neutral atoms 
in optical lattices [IF) and use the optimal control tech- 
niques to find efficient and robust pulse sequences for the 
logical couphng gate. 

We notice that other proposals which exploit the su- 
perlattice structure in 2D to create a universal resource 
(which is different from the universal cluster state) for 
MBQC by connecting Bell entangled pairs by entangling 
phase gate have been proposed recently |7]. 

The paper is organized as follows: first, in Sec. [IT] we 
present the general ideas for generating cluster states 
within a decoherence free subspace (DFS) using opti- 
cal superlattices. Then, in Sec. |III| we briefiy review the 
singlet DFS of the plaquette and describe operations of 
single logical qubit using superexchange couplings. Fi- 
nally, in Sec. |IV| we consider the key challenge of imple- 
menting the logical controUed-phase gate with the sin- 
glet DFS encoding. We propose three new approaches 
for the controUed-phase gate: the geometric phase ap- 
proach, the perturbative approach, and the optimal con- 
trol approach. A detailed comparison among the three 
approaches is summarized at the end of this section (Ta- 
ble [l| before we present our conclusion in Sec. [v[ 



II. DECOHERENCE FREE CLUSTER STATES 

One promising approach to quantum computation is 
the measurement-based quantum computation (MBQC) 
O [T5] , which uses universal resources such as the clus- 
ter states. The cluster states can be efficiently pre- 
pared by initializing all lattice spins in the product state 
of 1+) — (|0) -I- |1)) and performing the controlled- 
phase gate between all pairs of neighboring spins. The 
controUed-phase gate induces an additional factor —1 if 
both input qubits are in state |1). Up to some individual 
qubit rotations, the controUed-phase gate can be achieved 
by Ising-type interaction between two input qubits. 

The preparation of a cluster state has been demon- 
strated using optical lattices with logical qubits di- 
rectly stored in individual spins. ControUed-coUisions 
are used to implement the controUed-phase gate between 
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FIG. 1: (Color online) The optical superlattice consists of 
a periodic array of 2 x 2 plaquettes. The 2D optical trap- 
ping potential is created by adding two superlattice potentials 
Vx (x) and Vy (y) illustrated in the bottom and left panels. 
The intra-plaquette coupling is represented by the solid lines 
and the inter-plaquette coupling is represented by the dashed 
lines. 



neighboring spins. However, such atomic qubits for the 
controUed-collision scheme are vulnerable to magnetic 
field fluctuations, which limits the practical implemen- 
tation of the MBQC. 

The optical superlattice inducing a periodic array of 
2x2 plaquettes (as shown in Fig. [T|) can be created by 
superimposing two optical lattice potentials with short 
and long wavelengths differing by a factor of two [TB] 
along both, x and y, directions. The effective optical 
trapping potential becomes 



V^V, (x) + Vy (y) 



(1) 



where K = K.s cos^ (^u - (j)u,s) + K,; cos^ (f u - <j)u,i) 
for u — x,y. The short-lattice wavelength is A, and the 
parameters V^u,sj Vu,ii 4'u,s f^nd ; are controlled by the 
intensities and phases of the laser beams. 

For integer filling with one particle per site, each pla- 
quette has four particles. The four spin-1/2 particles 
have a two- fold singlet subspace with total spin zero along 
all directions (i.e., Stot — 0). Thus, the singlet subspace 
is the DFS insensitive to uniform magnetic field fluctua- 
tions. 

The intra-plaquette couplings (solid lines in Fig. [T]) en- 
able operations of single logical qubit encoded in the pla- 
quette. (Note that by manipulating intra-plaquette cou- 
plings minimum instances of topological matter can be 
demonstrated in the same optical superlattice |17j.l In 
order to create the decoherence free cluster states, we 
will also need inter-plaquette couplings (dashed lines in 
Fig. [T|) to implement the controUed-phase gates. 

There are eight sites for two neighboring plaque- 
ttes. If only four middle sites are involved for the 




FIG. 2: (Color online) Simultaneous coupling between neigh- 
boring plaquettes. (a) It is possible to simultaneously im- 
plement the controlled-phase gates between neighboring pla- 
quettes along the horizontal direction, if only four sites are 
involved for each controlled-phase gate (e.g., for the geomet- 
ric phase approach or the optimal control approach), (b) The 
optical superlattice (with inter-plaquette coupling and alter- 
nating energy off-set for odd and even plaquettes, created by 
A, 2A, 4A) can yield the inter-plaquette coupling for the geo- 
metric phase approach, (c) The spin-dependent optical super- 
lattice (indicated by thin and thick lines of potential profiles) 
can generate Ising interaction [13] useful for the optimal con- 
trol approach. 



controlled-phase gate (e.g., for the geometric phase ap- 
proach in section |IV A| or the optimal control approach 
in section 



IV C I , it is possible to simultaneously apply 



controlled-phase gates to couple all horizontal (or verti- 
cal) neighboring plaquettes with no overlap of sites in- 
volved for different controlled-phase gates as shown in 
Fig. [2] Meanwhile, if all eight sites from both plaquettes 
are involved for the controlled-phase gate (e.g., the per- 
turbative approach in section IV B I , two steps are needed 



to couple the plaquettes along the horizontal direction: 
first couple each even plaquette with the neighboring odd 
plaquette on the left, and then couple each even plaque- 
tte with the neighboring odd plaquette on the right. 

In order to use the prepared cluster state for the 
AIBQC, we should also be able to measure the individ- 
ual qubits. This can be achieved by first converting the 
spin singlet/triplet states into different particle number 
configurations [U [TB] , and then using various techniques 
of coherent optical control with subwavelength resolution 
[18} [T9| ^20^ to projectively count the particle number at 
a specific site (without compromising the coherence for 
the remaining sites). 

In the next two sections, we will consider the rotation 
of single logical qubit using intra-plaquette couplings, 
and the controlled-phase gate between two logical qubits 
using the additional inter-plaquette couplings, respec- 
tively. 
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of total spin 0. We use the two-fold singlet subspace of 
the plaquette to encode the logical qubit. The singlet 
subspace is spanned by 



\^h) = \S),,(^\S) 



3,4 



(4) 




FIG. 3: (Color online) Intra-plaquette superexchange cou- 
plings. The coupling strengths are (a) Jh = t'n/U between 
horizontal neighbors and (b) Jv = iy/U between vertical 
neighbors. They can be changed independently by tuning the 
barriers between the sites. 



III. LOGICAL QUBIT ENCODED IN THE 
PLAQUETTE 

In this section, we focus on the operations within the 
plaquette via intra-plaquette coupling. For concreteness, 
we consider bosonic particles, and similar results can be 
obtained for fermionic particles as well. The following 
Hubbard Haniiltonian governs the dynamics of a single 
plaquette, with spin-independent tunnelings and interac- 
tions (we will introduce an additional vibration level in 
Sec. IV): 

H = - ^ {tija\^aj„ + h.c^ + ^[/^ n, (ni - 1) 
+ '^^J,ini^ (2) 

i.(T 

where (al^) is the annihilation (creation) operator, 
rtiCT is the particle number operator for site i = 1, • • • ,4 
with spin a ='\,ij ^^J^d rii = n^.i -|- n.^^j^. The tunnel- 
ing amplitudes (i/f = ti2 — t^^ and ty = = ^41) 
and the offset energies {fj,i) can be changed by tuning 
the superlattice parameters. The large on-site interac- 
tion U ^ tH,tv ensures that the system is in the Mott 
insulator regime with fixed particle number for each site. 

Particle tunneling only occurs virtually between neigh- 
boring sites, which leads to the superexchange coupling 

Heff = - Jh (Sl ■ 5*2 + S3 • S4) - Jv {S2 • S3 -I- S4 • Si) , 

(3) 

where si are Pauli operators for the spin at site i. The 
coupling strengths Jh = tn/U and Jy = ty/U can be 
changed independently, by tuning the barriers between 
the sites as illustrated in Fig.|3] Controlling the superex- 
change couplings is sufficient to perform arbitrary rota- 
tions of the logical qubit encoded in the plaquette. 



A. Singlet subspace for four spins 

The space of 4 1/2-spin particles span a subspace of 
total spin 2, 3 subspaces of total spin 1 and 2 subspaces 



\^v) = \S),,^\S) 



/4,1 ' 



(5) 



with \S)^^^ ^ ^ (it), U),. - li), IT),). \^h) (or \^v)) is 
the product state of two singlet pairs along the horizontal 
(or vertical) direction, which can be prepared using the 
procedure demonstrated in 

The singlet subspace is decoherence free, because it 
is insensitive to the uniform magnetic field fluctuations. 
In addition, measuring a single spin will not distinguish 
the states from singlet subspace, and this is a source of 
protection against local perturbations. Since (^'^fl^'v) = 
1/2 0, it is more convenient to use the orthogonal 



states |0) = and |1) 
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V3 



l^'y) - I^'h))- We can 



also write the orthogonal states in terms of the singlets 
and triplets for vertical pairs (2, 3) and (4, 1) pT] : 



1^)2,3 l'^). 



(6) 
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|7^+)2,3l^^-)4,l 



^0)2,3 1^0) 
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(7) 

where |r+,o.-) - {|T) IT) , (IT) U) + li) IT)) , li) li)}- 
For such a choice of basis, the subsystem of two spins 
(2, 3) is sufficient to determine the logical states |0) and 
|1), because the corresponding reduced density matrices 



i\S) (5|)2,3 
1 



Tr4,i [|0) (C 

Tr4,i[|i)(i|]-^(|r+)(r+| + |ro) (To\ + \ 



(8) 
(9) 



belong to orthogonal (singlet and triplet) subspaces. 



Tr 



^'2,3/^2, 3 



0. 



The Pauli operators associated with the logical qubit 
are: a, = |0) (1| |1) (0|, ay = i\0) {1\ - i\l) (0|, and 
az = |0) (0| — |1) (1|. Within the singlet subspace the 
logical operator can be achieved by operating the (2, 3) 
spins 



1 



(1 + S2 • S3) . 



(10) 



where we use " =" to represent the special equality valid 
within the singlet subspace. Since 32 • S3 = S4 • si (i.e., 
either both pairs are singlets or both are triplets) , can 
also be implemented by operating the (4, 1) spins 



1 



(H- S4 • Si) . 



(11) 
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For example, |-| 
two-step evolution e 



FIG. 4: (Color online) The Bloch sphere representation for 
the singlet subspace. The state |0) (|1)) is associated with the 
north (south) pole, (a) The superexchange coupling Jh (Jv) 
is associated with the rotation around the uh (nv) axis. Here 
uh ~ (^^lOj^^ S'lid fiv = (0,0,1). The sequential rota- 
tions around the uh and nv axes can rotate the Bloch vector 
from fiv to rig = {^,0, ^) [i.e., state ^ (|0) -I- |1))]. (b) 
Alternatively combined superexchange coupling (with contri- 
butions from both Jh and Jv) can implement the rotation 
around the axis nc, which rotates the Bloch vector from nv 
to fig in one step. 



B. Rotating logical qubit with super-exchange 
coupling 

We now consider arbitrary rotations in the sin- 
glet subspace using super-exchange couplings. First 
of all, the superexchange coupling Hamiltonian com- 
mutes with the total spin operator of the plaque- 
tte [Hf,ff, si -I- 5*2 + S3 -t- S4] = 0, due to the identity 
[si ■ Sj , Si + Sj] = 0. Consequently, the superexchange 
coupling preserves the singlet subspace (with zero total 
spin). Within the singlet subspace we have 

n/^ • (T ^ i - i (si • S2 + S3 • S4) (12) 

• (T = i - i (52 • S3 + S4 • Si) (13) 



with fiH — (^^,0, and ny = (0,0,1) as illustrated 

in Fig. Qa). The constant of 1/2 can be neglected, as it 
only induces an overall phase during the evolution. The 
rotations about these axes can be controlled by switching 
on/off the superexchange couplings of Jh and Jy, which 
varies expentially with the height of the corresponding 
barriers. Since the angle between nn and ny is 27r/3, 
arbitrary rotation of the Bloch sphere can be achieved 
within 4 operations. (This is a special case of the general 
theorem [5TJ|52] stating that k+2 operations are sufficient 
for arbitrary rotation given the angle rj between the two 
rotation axes satisfies ^ > min (77, tt — 77) > j^-) 

The product state of two vertical singlet pairs |0) — 
l^'y) can be initialized using the procedure demonstrated 
in [3]. Universal rotation enables dynamical preparation 
of arbitrary logical state encoded in singlet subspace. 



2 sin 



1)) can be prepared by 
^-y) with 20H = 
sm - ^ « 125.3°, as 



^ « 109.5° and 26'y = Tr-sin"^ ^ 
shown in Fig. [IJ^a). Alternatively, we can tune the rela- 
tive strength between n/f •(? and ny •(? to achieve the total 

coupling nc'CT with — and prepare |+) in 

one step e"*"'^''^? [^-y), as illustrated in Fig.|4];b). Note 
that all the plaquettes can be simultaneously prepared 
in the |-|-) state. In order to create the decoherence free 
cluster state, we need the controUed-phase gate between 
the logical qubits encoded in neighboring plaquettes. 



IV. CONTROLLED-PHASE GATE 

We now consider inter-plaquette couplings (dashed 
lines in Fig. [ij. In particular, we focus on implement- 
ing the controUed-phase gate between two neighboring 
plaquettes, which induces an additional —1 phase if both 
encoded qubits are in the logical state |1). In princi- 
ple, the controlled-phase gate can be achieved by the 
Ising-type interaction between the logical qubits, but un- 
fortunately such interaction is not immediately available 
from the lattice experiments, as the effective Ising term 
(JzO'z requires four-site interaction (5*2 • S3) (s^ • s'^) [see 
Eqs. ( 10|11 1]. However, since what we want is the spe- 
cific unitary evolution rather than the interaction, it is 
actually more feasible to implement the unitary evolution 
directly. 

In the following, we present three different approaches 
to implement the controlled-phase gate between two 
neighboring plaquettes. For concreteness, we only con- 
sider coupling two neighboring plaquettes along the hori- 
zontal direction, while all three approaches can also cou- 
ple neighboring plaquettes along the vertical direction. 



A. Geometric phase approach 

The first approach uses the vibration levels and the 
geometric phase to achieve the controlled-phase gate be- 
tween neighboring plaquettes [Fig. [5]ja)]. The geometric 
phase is proportional the surface area enclosed by the 
evolution trajectory in the Bloch sphere (associated with 
the two energy levels that are degenerate). For example, 
if a half of the Bloch sphere is enclosed, the system ac- 
quires a geometric phase tt. We first consider the bosonic 
particles. It takes three steps to achieve the controlled- 
phase gate: 

Step 1. We lower the inter-plaquette barrier and adia- 
batically tilt the intra-plaquette potential along the verti- 
cal direction [Fig. [sjb)]. Each lower site will be occupied 
by one particle (or two particles) if the vertical pair of 
particles is in the singlet (or triplet) state [Fig. [5][c)]. For 
example, if the spins (2,3) are in the singlet state (de- 
noted as 52,3), the transfer of particle from site 2 to site 
3 is prevented by the symmetry requirement of bosonic 
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FIG. 5: (Color online) Geometric phase approach to 
controlled-phase gate (Step 1). (a) The sites from two neigh- 
boring plaquettes are labeled, (b) The intra-plaquette trap- 
ping potential along the vertical direction is adiabatically 
tilted. This results in single (or double) occupancy at the 
lower site if the vertical pair of particles is in the singlet (or 
triplet) state, (c) Particle number configurations are plotted 
for four possible of spin states: ^2,3 ® Si/_4/, 82,3 <8> Ti/,4/, 
T2.3 <8) 5*1/, 4', and T2^3 (8) Ti/_4/. (S and T indicate singlet and 
triplet.) 



particles, resulting in one particle in site 3 [see the up- 
per two panels in Fig. [sjc)]. If the spins (2, 3) are in the 
triplet subspace (denoted as 12^3), the particle from site 2 
is adiabatically transferred to site 3, leaving two particles 
in site 3 [see the lower two panels in Fig. |5jc)]. Similar 
spin-dependent transfer also happens to other sites, such 
as (l',4'). 

Step 2. We quickly apply a defined bias to the inter- 
plaquette lattice potential and lower the inter-plaquette 
barrier along the horizontal direction [Fig. [6]^b)]. This 
induces single particle resonant tunneling with rate t be- 
tween the vibrational ground state at site 2 and the vibra- 
tional excited state at site 1' [33 , if there is one particle 
at site 2 and zero particle at site 1' [see the highlighted 
upper right panel in Fig.|6jc)]. By waiting for time 27r/i, 
we obtain the geometric phase tt from the resonant tun- 
neling for S'2,3<8)ri/.4/. As detailed in Sec. IV, for all other 
three cases (6*2,3 ® •S'l/^', T2,3 (g) S'1',4', and T2,3 (g) T'1/,4/) 
we only obtain a trivial geometric phase (or 2tt) . 

Step 3. We change the intra-plaquette potential to 
the initial balanced position along the vertical direction 
(having one particle per site) and restore each plaquette 
to the logical subspace. 

A recent superlattice experiment uses the resonant 
tunneling and the blockade induced by on-site interac- 
tion to count the number of atoms |23j . This experiment 
demonstrates that the presence or absence of resonant 
tunneling can be highly sensitive to the number of par- 
ticles in the lattice sites. The geometric phase approach 
can be regarded as an extension that uses the resonant 
tunneling to coherently imprint a geometric phase for a 
specific particle number configuration (corresponding to 
certain logical state). 

The procedure for the fermionic particles is almost 
the same as that for the bosonic particles, except for 




FIG. 6: (Color online) Geometric phase approach to 
controlled-phase gate (Step 2). (a) The sites from two neigh- 
boring plaquettes are labeled, (b) A defined bias (A) of 
the inter-plaquette potential is quickly applied and the inter- 
plaquette barrier is lowered to facilitate the resonant tunnel- 
ing along the horizontal direction, (c) Resonant tunneling 
(between the vibrational ground level of the left site and the 
vibrational excited level of the right site) can occur for the 
following two cases: 1) each of the left and right sites has 
exactly one particle, and the two particles are in the singlet 
state (see the left panel), 2) the left site has one particle and 
the right site has zero particle (see the highlighted upper right 
panel). All other configurations are off-resonant, with negli- 
gible tunneling. After time 27r/t, a geometric phase tt from 
the resonant tunneling is obtained for 52,3 (8 ri',4' (the high- 
lighted upper right panel) , while only a trivial geometric phase 
(0 or 2-k) is obtain for the other three cases (52,3 (8) Si/. 4/, 
12,3 (8 S'1',4', and T2,3 (g)ri',4/). 



the following three differences. First, the bias of the 
energy off-set needs to be A = + If^ for fermionic 
particles (whereas A = lo for bosonic particles), where 
uj is the vibrational excitation energy and is the 
on-site interaction between ground and excited levels at 
the right site (1' or 4'). Second, the geometric phase tt 
is obtained from the resonant tunneling associated the 
subspace T2,3 S'1',4' for fermionic particles (whereas 
it is associated with ^2,3 T'i',4' for bosonic particles). 
Third, the geometric phase is for the remaining cases 
for fermionic particles (whereas it might be either or 
27r for bosonic particles). 

It is tempting to consider using A = u for the fermionic 
particles, as we might expect that by exchanging the roles 
of singlet and triplets, the fermionic particles could be 
mapped to bosonic particles. However, the roles of singlet 
and triplets are not exactly symmetric. For example, 
consider the case with one particle per site after Step 
1. For bosonic particles, the system is in the subspace 
62,3 (X" S'1',4' that has finite projection to S2,i' (E" 63,4' and 
^2,1' (E> T3,4' (but not r2,i' (g) S3,4' or 62,1' 13,4'), which 
yields a trivial or 27r geometric phase. For fermionic 
particles, the system is in the subspace r2,3 ® Ty^' that 
has finite projection to T2,i' ® 63,4' and S2,i' ® T^a' , as 
well as S2,i' (8) 63,4' and T2.1' (g) T3,4', which thus may 
yields a non-trivial tt geometric phase. 

The detailed calculation for both bosonic and fermionic 
particles are presented in Appendix \K\ 
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B. Perturbative approach 

The second approach uses both the intra- and inter- 
plaquette couplings acting on the eight sites. The intra- 
plaquette coupHng induces an energy gap between the 
logical states (i.e., singlet subspace) and other non-logical 
states, while the inter-plaquette coupling acts as a per- 
turbation that induces different phase shifts for differ- 
ent logical states. The inter-plaquette coupling can be 
efficiently achieved using super-exchange interaction be- 
tween the inter-plaquette neighboring sites. 

The key challenge is to obtain the intra-plaquette in- 
teraction, with finite Heisenberg interaction between the 
sites along the diagonal and off-diagonal directions. We 
can overcome the challenge by using a different design of 
the optical lattice. 



1. Lattice Geometry and Energy levels 

We want to obtain the Hamiltonian with intra- 
plaquette interaction: 



Hi, 



= J 



E 

=1,2,3,4 



Si ■ Si-l-2, 



(14) 



where the exchange interaction J — ±.t^ /U {d — zLP /U) 
is induced by the tunneling between the nearest neigh- 
bors (next-nearest neighbors) with tunneling rate t (t). 
The positive and negative signs are for fermions and 
bosons, respectively. To simplify the notation, we have 
identified S5 with si. Such type of interaction can be 
created by a lattice potential of the form [sec Fig. [7][a)] 



□ 



A 4 - 




» 





=0^ S=1^ S=3| 


4{ZJ+d] 


|2> 




, , |1H,«1> 


^ 4lJ-d| 




4d 


4 ► M'"> 


-4(J-d) 


— !□> 



FIG. 7: ( Col or online) (a) Density plot of the lattice potential 
[see Eq. ( 15 1], which generates an array of 2 x 2 plaquettes in 
the X — y plane with periodicity of 2n/k. Using parameters 
= eO^:^, V = 9Sr and V' = IQEr, with Er and Er pho- 
ton recoil energy of the long and short lattices respectively, 
one can achieve a parameter regime with d/J ~ 0.2. (b) En- 
ergy levels of a single plaquette described by the Hamiltonian 
given in Eq. (14 1. In the plot we assume fermionic atoms, i.e. 

For bosons J, d < the order of the en- 
2 is the lowest energy state. 



J,d > (fermions) 
ergy levels is reverse i.e the 5* - 



where V^{u) = ^cos(2fcw) — ^cos(fcu) are the typi- 
cal double well superlattice formed by the superposition 
of two independent sinusoidal potentials which differ in 
periodicity, a = 7r/fc, 27r/fc, by a factor of two 3 and 
V'^{x,y) — —V'^{k{x — y))cos'^{k{x + y)) is an additional 
potential that allows to control the diagonal couplings 
within the plaquettes. ft can be constructed for example 
from a folded, retro-reflected beam with out-of-the-plane 
polarization [24]. By varying the depths of the short, V 
and long lattices it is possible to control the intra and 
inter-plaquette coupling independently, and in particular 
to make the latter negligibly small and the plaquettes 
independent. 

As the intensity of the non-separable part of the poten- 
tial is ramped up, a minima at the center of the plaque- 
ttes develops. If the strength of the latter is such that 
the energies of bound states in this minima are larger 
than the energies of the lowest vibrational states at the 
plaquette sites it is possible to tune the ratio t/t with- 
out populating the central site, which is required for the 
validity of Eq. ( |14[ ). For example using the parameters 
= 60Er = 15Er, V = 9Er and V = lOEr, with 
En = fi^k^/iSm) and Er = H'^k'^/{2m) the photon re- 
coil energy of the long and short lattices respectively one 
can achieve a parameter regime with t/t « 0.5 with an 
energy gap to the first vibrational state in the central 
well of order Eg /t Ri 10. ft is very difficult to increase 
t/t close to one by just controlling the lattice potential, 
because the energy gap disappears and the central sites 
become accessible. Therefore, we will focus on the case 
d<J. 



The eigenstates associated with Eq. (14 1 can be clas- 
sified according with their total spin S. As shown in 
Fig.|7jb), there are two singlet {S = 0) states 



\a)^^{\^H) + \^v)) 

|X> = - l^-y). 



(16) 
(17) 



V{x, y, z) = V\x, y) -h V^{x) + V^{y), (15) 



with energies E{\U)) = -4{J - d) and £;(|x)) = +4( J - 
d), respectively. There are three S = 1 states denoted 
by |l(-i'0'i)), with energies E{l'^i^) = 4J,4J and Ad, re- 
spectively. For fermionic (bosonic) atoms the highest 
(lowest) energy state is a 5' = 2 state, |2) with energy 
E{2) = A{2J + d). 

We want to use the singlet states within each pla- 
quette as encoded qubits and perform a phase gate be- 
tween them by coupling nearest neighbor plaquettes into 
a "superplaquette^' (i.e., 2x4 potential wells). A super- 
plaquette can be achieve by superimposing laser beams 
with periodicities 4 A and 4 A/3 along one axis [25 . Such 
wavelength are experimentally available for typical Alkali 
atoms or can be engineered by intersecting pairs of laser 
beams at appropriated angles [55]. The 4A isolates pairs 
of adjacent plaquettes along one direction and the extra 
4 A/3 lattice is needed to balance the offset created when 
the latter lattice is added. When pairs of plaquettes are 
weakly coupled into a superplaquette the Hamiltonian 



that connects the plaquettes is given by 

He = J'{S2 ■ Sy + S3 • S4/). 



(18) 



We want to use the coupHng to implement a controUed- 
phase gate between the singlet eigenstates in the 
two plaquettes. To achieve that we require that 
the inter-plaquette coupHng is weak (i.e., J' <C 
min{4c?, 8 ( J — d) , 4( J — 2c?)}) and derive an effective 
Hamiltonian by adiabatically eliminating the all S* > 
states. 

In the following we discuss the implementation of 
the controUed-phase gate for the experimentally relevant 
regime d < J. (The ideal case of d = J is discussed in 
Appendix Ib] ) 



2. Perturbative approach with d < J 

For the 2D plaquette implementation the diagonal cou- 
pling d is always smaller than J and therefore the singlet 
states within the plaquette are non-degenerate 



|A£;| = |i?(|x))-i?(P))| = 8|J-d|>0. 



(19) 



Regardless of this issue, it is still possible to derive an 
effective Hamiltonian provided that the inter-plaquette 
coupling J' is less than the energy difference between 
|x) and all other states [see Fig. [jf^b)]: 



J' < min{4d, 8 ( J - d) , 4( J - 2d)}. 



(20) 



From this consideration we observe that close to d = 
0, J/2 and J, the perturbative approach (based on !□) 
and |x)) breaks down and we should stay away from 
these points. 

As detailed in Appendix [B] we can obtain the effective 
Hamiltonian 



zreff _ , 



J 



) E -I 



(21) 



•11 
J 



-ffi ■ Si+i + (A^ - -)o-f (jf+i 



where a are effective Pauli matrices acting on the !□) , | x ) 
states, and 



9J 




8J 


24J 


T ^ 


d 


-3J 


+ 2-3 + 

d+ J 


9J 




8J 


^ 2J-d) 


T + 


d 


-3J 





J 



(22) 
(23) 



In Fig. |8]^a) the parameters A^ and 7^ are plotted as a 
function of d/ J. 

Within a superplaquette the term Ui ■ (Xi+i commutes 

with Hl{^-^ (it only introduces a phase 4>t = ^jg^c in the 
effective triplet subspace: ]□,□) ,|x,x) and ([□, x) + 
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FIG. 8: (Color online) (a) Parameters of the effective Hamil- 
tonian for the general case d ^ J. At the point d ~ 0.62J , 
the Ising term in the effective Hamiltonian vanishes, = | 
(see text), (b) The controUed-phase gate fidelity F on a su- 
perplaquette as a function of d/J for different rations of J'/ J. 
For same d/J, the smaller J'/J the higher the fidelity. There 
are 4 critical points at which the fidelity drops considerably, 
these are d/J « 0,0.5,0.62,1. At d/J « and 0.5, one of 
the singlet states becomes degenerate with one S — 1 state 
and consequently the effective Hamiltonian breaks down. At 
the d/J « 0.62 the Ising term vanishes and at d = J the 
rotating wave approximation used in the simplification of the 
effective Hamiltonian becomes invalid. The shadow regions 
are the ones where the achievable fidelity is higher than 0.98. 



\x,D))/y/2 and (/)5 = "ijjs^c for the effective singlet 
states: (]□, x) — |x,n))/\/2). Here tc stands for the du- 
ration of the controUed-phase gate, i.e., ^tc{\z ~ \) — 



h{2n- 



h: 



eff 



l)7r/4, with an integer n = 1, 2, • • • . Consequently 
can be used to perform a controlled-phase gate 
within a superplaquette. We use the standard echo tech- 
nique (i.e., TT pulses at tc/2 and tc for each of the encoded 
qubits) to remove the unwanted ctJ term from the effec- 
tive Hamiltonian in Eq. (21 1. The controlled-phase gate 



can be achieved by the unitary evolution 

U ^ X e-i("i^t--+Hc)tc/-2 X (,-iiH,„t,-a+H,)tj2^ ^24) 

where X represents the echo tt pules for the encoded 
qubits, which can be achieved via intra-plaquette super- 
exchange couplings. 
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FIG. 9: (Color online) For the d < J, the cluster state gener- 
ation has to be applied in two steps (a) and (b). 



We use the exact diagonalization to calculate the 
controUed-phase gate fidelity F — |/|^ with 



/ = — Tr 



U. 



c— phase 



PUP 



(25) 



where P is the projection operator to the singlet subspace 
for each plaquette. In Fig.jsjjb), we plot as a function of 
d/J for different rations of J' /J, with n — 1. The Figure 
shows that at the points d w 0.5J, d « 0.62J, d = 0, 
d = J there is an abrupt drop of the fidelity, as the ratio 
J' /J is increased. The drop at these points is expected 
since at d « 0.5J and d « one of the singlet states 
becomes degenerate with one S = I and consequently 
the effective Hamiltonian breaks down. At c? « 0.62J 
, Az = 1/8 [see Fig. [sj^a)] the Ising term vanishes and 
at d = J the two singlets become degenerate and the 
rotating wave approximation (assumed for the derivation 
of the effective hamiltonian in Appendix |b]) is not longer 
justified. 

Away from these points the derived effective Hamil- 
tonian provides a good description of the dynamics and 
for values of J' w O.IJ one can get a gate fidelity above 
0.98. In the plot we highlight with a gray shadow the 
d/J parameter regime where the achievable fidelity is 
above 0.98. However among these shadow regions only 
the regime d/J < 0.5 is experimentally achievable us- 
ing the lattice geometry described early in this section. 
The small fiuctuations in the fidelity curves are due to 
the non-energy preserving terms neglected to obtain the 
effective hamiltonian (see Appendix |B]), which can be 
suppressed when J' <C ^^8 (J — d) J. 

The Heisenberg term CTj • (Ti+i, on the other hand, does 
not commute with H^^( consequently the phase 

gate can not be applied simultaneously to all plaquettes. 
Instead it has to be applied first between the superpla- 
quettes formed by the plaquettes 2i + 1 , 2i + 2 and sub- 
sequently between the superplaquettes formed from pla- 
quettes 2i + 2,2i + 3 (see Fig. [9]). Additionally in order 
to create a cluster state across all the plaquette array, it 
is required to fine tune the parameters and time evolu- 
tion to eliminate the different phase accumulated by the 
triplet and singlet states in the encoded spin basis due to 
the Heisenberg term at tc- Consequently for multipartite 



entanglement generation not all d/J values are allowed 
but only the ones which satisfy the following conditions: 



(pT — 4>s — 27rm 



^tc{K-\) = h{2n-l)'^, 



(26) 
(27) 



where n and m are integers. In Fig. 10 'a), we show a 
set of allowed d/J values which satisfy the conditions 
given by Eq. (26 1, for different n and m values. Here 



we also highlight with a gray shadow the correspond- 
ing d/J values which yield a fidelity higher than 0.98 for 
J' / J < 0.1. In Fig. 10 b) we show two examples of traces 
of the phase gate fidelity vs J' / J:{n,ni) = (1,1), (3, 4) 
(indicated in panel a. by a square) computed by exact 
diagonalization of the superplaquette Hamiltonian. The 
figure shows that it is always possible to find parameters 
which allow for a high gate fidelity. However, here we 
are only including errors due to higher order terms ne- 
glected in the derivation of the effective Hamiltonian. In 
realistic experiments other external errors such as lattice 
inhomogeneities are always present, which can be mini- 
mized at the expense of a larger J' /J ratio (faster evo- 
lution) . There is consequently a trade off between faster 
time evolution and small perturbative corrections. 

In conclusion we have presented a scheme to perform 
controUed-phase gates in the encoded singlet subspace. 
This perturbative scheme has two advantages: (1) the 
plaquettes are always in the decoherence free subspace, 
(2) it is easy to implement as it only relays on the coher- 
ent dynamical evolution without further manipulations. 
Due to the fact that the dynamics is determined by a 
second order effective Hamiltonian the achievable fidelity 
with the proposed schemes can become very high but at 
the cost of slower time evolution. If the strongly interact- 
ing regime is reached by using a Feshbach resonance, one 
can achieve values of J of order of 100 Hz and therefore 
cluster generation times of order 0.1 — 1 sec. These gener- 
ation times are slow but longer than the encoded qubits 
decoherence time due to their insensitivity against envi- 
ronmental decoherence [3]. 



C. Optimal control approach 

We now consider the optimal control approach to fast, 
high fidelity implementation of the controUed-phase gate 
[between the horizontal neighboring plaquettes (1,2,3,4) 
and (l',2',3',4') as shown in Fig. ^a.)]. The key chal- 
lenge here is to identify an efficient set of operators that 
(i) enable the unitary evolution of the controUed-phase 
gate, and (ii) are feasible using optical superlattices as 
well. We first provide a set of operators sufficient to 
achieve the controUed-phase gate with arbitrary preci- 
sion. After that we numerically find the pulse sequences 
for these operators to implement the controUed-phase 
gate. 
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FIG. 10: (Color online) Examples of controlled-phase gate 
with perturbative approach, (a) The values of d/J that sat- 
isfy the conditions stated in Eqs. (26 1. The shadow re- 



gions highlight the regime where the phase gate fidelity can 
be larger than 0.98 for J'/J < 0.1. (b) The controlled-phase 
gate fidelity f as a function of J'/J can be calculated by ex- 
act numerical diagonalization. The red dashed (or blue solid) 
curve is for (n, m) — (1,1) (or (3,4)), which is also indicated 
by the red (blue) square in panel (a). 



Choice of operators 



for the set of operators 

C2 = • S4 



O3 = S2,zSi,2 + 53,^54,2 



(29) 
(30) 
(31) 
(32) 
(33) 



To justify that H (t) can implement the controlled- 
phase gate 



c-phasc — GXp 



-z-(l-a,)(l-a;) 



(34) 



we show that ( 1 — ct^ ) ( 1 — cr^, ) belongs to the Lie alge- 
bra generated by {Ofej^^i ... 5- We start with these five 
operators as the available set (AS), and calculate the 
commutators among the AS operators. We then expand 
the AS by adding new commutators that are not lin- 
ear combinations of the AS operators. We denote the 
number of linearly independent AS operators as the di- 
mension of the AS. We repeat the process of calculating 
the commutators and expanding the AS, until its dimen- 
sion does not increase any more. We use Mathematica 
to iterate the process of expanding the AS until it sat- 
urates at dimension 80 (including the identity operator 
that commutes with all other operators). Finally, we 
verify that (1 — cr^) (1 — a'^) (5*2 • S3) (s\ • s^) is a lin- 
ear combination of the AS operators. Therefore, accord- 
ing to the local properties from the Lie algebra the set 
of operators {Ofc}j.=i ... 5 is sufficient to implement the 
controlled-phase gate. The remaining task is to find the 
solution for {ak (t)}- 



In principle, tunable Heisenberg superexchange inter- 
actions are sufficient for the controlled-phase gate by cou- 
pling all eight sites pH. However, for the optical lat- 
tice experiments, we would like to achieve the controlled- 
phase gate by coupling as few sites as possible, preferably 
using global rotations for all spins and Heisenberg/Ising 
interactions between neighboring sites |5] |T3] . 

According to Eqs. ( 10|11 ), we need at least two sites 
from each plaquette to determine the cr^ operator. Since 
the controlled-phase gate depends on both operators 
from the plaquettes, we should consider at least four sites 
to implement the controlled-phase gate. It turns out that 
coupling the four middle sites (2,3,4', 1') is sufficient to 
achieve the controlled-phase gate, which significantly re- 
duces the complexity compared with the earlier proposal 
that couples all eight sites [TT| . 

We consider the Hamiltonian 



iI(t) = ^«fe(t)Ofc, (28) 
where {ak (t)} are the time-dependent control variables 



2. Smooth pulses 

We use an algorithm which can be interpreted as a con- 
tinuous version of the gradient ascent pulse engineering 
(GRAPE) [271 US], though it is developed in an indepen- 
dent way |29j . The algorithm based on optimal quantum 
control is summarized in Appendix [C] Comparing with 
the GRAPE method it has the advantage that we can find 
solutions with specific boundary conditions (e.g. pulses 
start and end at zero) and in terms of smooth (finite 
slope) functions of time as well. 

More specifically, the particular form of the coefficients 
{ak (t)} are chosen to be finite sums of sinusoidal func- 
tions 



ak{t, Xki, ■ . ■,XkL) = '^Xki sin 
1=1 



lTTt\ 



(35) 



each of which depends on L parameters {xki\ (I — 
1,...,L). As mentioned before, note that they fulfill 
the convenient property that Q!fc(0) = ak{T) = and 
that they have a finite slope. Hence, we need to optimize 
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Geometric Phase 


Perturbative 


Optimal Control 




Approach 


Approach 


Approach 


Time scale for controUed-Z (CZ) gate 


1/J 




1/J 


Duration out of DFS 


lit 





1/J 


Systematic errors 


(t/Uf 


{J'/Jf 





Inhomogeneity errors 


{5t/tf ,{5A/tf 


{5J/Jf ,{5J'/J'f 


{5J/Jf 


Sites per CZ gate 


4 


8 


4 


Simultaneous coupling 


Yes 


Two steps 


Yes 


A/TflirtT" iTitprflptimm 


Superexchange 
Single particle tunneling 


ouperexchange 


Superexchange 
Ising interaction 


Superlattice wavelengths 


A,2A,4A 


A,2A,4A,4A/3 


A,2A 


Vibrational levels 


Ground + excited 


Ground 


Ground 


Physical process 


Clear 


Clear 


Hard to interpret 


Control complexity 


Medium 


Low 


High 



TABLE I: Comparison among three approaches. See Sec. |IV Dj for discussions. 



K X L (with K = 5 operators in our case) parameters 
Xki which maximize the fidelity F = where 



/=^Tr[C/lp,_C/(T;a:)] 

N 



(36) 



for a particular subspace of N states {IV'n)} ("■ = 
1, . . . , N) of dimension d> N. 

In Fig. 11 'a) we show an example of pulses obtained 
with i = 20, which can attain very low infidelity e — 
1 ~ F, less than 10"'' (the value can be further reduced 
by improving the precision of the numerics). In Fig.[Tl|b) 




FIG. 11: (Color online) (a) Smooth pulses with L — 20 yield- 
ing to infidelities smaller than 10^^. (b) Infidelity e = 1 — F 
v.s. the deviation 5 J/ J for the smooth pulses. The infidelity 
is constant at 10~^ for 5J/J < 10~* and increases quadrati- 
cally with SJ/J for 5J/J > 10"''. (c) Quadratic dependence 
of the fidelity with SJ/J. 



we plot the infidelity e as a function of the relative devia- 
tion SJ/J. We assume that the couplings in Eq. (28) are 
deviated from ak to (l — SJ/J)ak', that is the system 
evolves under the deviated Hamiltonian (1 — &J/J)H. 
For simplicity, we consider the case that 5J / J is time- 
independent. (For example, imperfect calibration of bar- 
rier height or barrier thickness may induce such pro- 
portional, time-independent deviation in superexchange 
couplings.) We find that the infidelity remains a con- 
stant value (approximately 10~^) for very small devi- 
ations (with &J/J < 10"''), while the infidelity scales 
as {SJ/J) for larger deviations (with SJ/J > 10"^) 



which is also plotted using the linear scale in Fig. 11 'c). 
Such quadratic dependence to the deviation is not un- 
common, as the infidelity for single spin rotations also 
scales quadratically with the deviation. The quadratic 
dependence can be regarded as a direct consequence of 
the optimization procedure, which finds a local minimum 
of the function with first order derivatives being zero. 



3. Experimental implementation 

We now briefiy discuss the implementation of opera- 
tors {Ok}k^i ... 5 (and {Ok}) for the cluster state prepa- 
ration. The operators of Oi and O2 can be achieved by 
super-exchange interaction using superlattice techniques 
[3], while the operator of O3 can be obtained from spin- 
dependent tunneling in optical lattices |13j . Further- 
more, we note that the evolution of the Ising interac- 
tions O3 between all horizontal neighboring plaquettes 
can be performed simultaneously, because they act on 
different groups of physical spins as illustrated in Fig. |2] 
For the same reason, the operators of O4 and O5 can 
be performed simultaneously for all spins by driving the 
entire optical lattice with appropriate micro-wave pulses 
[3]. Therefore, the simultaneous controUed-phase gates 
between all horizontally neighboring plaquettes can be 
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achieved. 



CONCLUSION 



D. Comparing three approaches 

We compare the three approaches (see Table |l]) in the 
following aspects: (1) two relevant time scales: the time 
to implement the controUed-phase gate, and the duration 
for the plaquette not being protected by the DFS (which 
should be short compared to the coherence time outside 
the DFS [34), (2) two types of errors contributing to 
the controUed-Z (CZ) gate infidelity: the systematic er- 
rors from the approximations used in our analysis, and 
the inhomogeneity errors due to the fact that the cou- 
plings (e.g., t, A, J, and J') are not exactly the same 
for all plaquettes, (3) the number of sites involved for 
each controUed-phase gate: if each gate only couples 4 
sites, the controUed-phase gates between all horizontal 
(or vertical) neighbors can be achieved simultaneously; 
otherwise two sequential steps are needed, (4) the ma- 
jor interactions, such as singlet particle tunneling |16j . 
(Heisenberg) superexchange coupling [3^, and Ising in- 
teraction |13j , (5) the wavelength components needed to 
construct the super lattices, (6) the relevant vibrational 
levels, (7) the interpretation of the physical process, and 
(8) the control complexity for time-dependent parame- 
ters. 

The maximum achievable fidelity for the CZ gate 
is limited by the systematic and inhomogeneity errors. 
For the geometric phase approach, the systematic error 
{t/U)^ is due to the off-resonant tunneling, which is ana- 
lyzed in Table |TT] and Table |III| for bosonic and fermonic 
particles, respectively; such off-resonant tunneling can be 
suppressed by using Feshbach resonances to increase U 
while keeping the same tunneling rate t. For the per- 
turbative approach, the controUed-phase gate fidelity is 
F > 0.98 for J'/ J < 0.1 (Fig.fTol. For the optimal con- 



trol approach, the systematic error is only limited by the 
precision of the numerics [Fig. 11 'b,c)]. We note that it 
is important to suppress the inhomogeneity errors, as all 
three approaches are sensitive to such imperfections. Re- 
placing the parabolic trap with the flat-bottom trap [30] 
can be one possible solution to reduce the inhomogene- 
ity errors. It would also be interesting to consider other 
approaches that are insensitive to the inhomogeneity er- 
rors. 

Overall, the geometric phase approach has the advan- 
tage of fast operational time, short unprotected duration, 
and compatibility of simultaneous coupling. The per- 
turbative approach has the advantage of always being 
protected by the DFS and favorable control complexity. 
The optimal control approach has the advantage of fast 
operational time, vanishing systematic errors, and com- 
patibility of simultaneous coupling. 



In conclusion, we have discussed preparation of large 
cluster states for neutral atoms in optical super lattices. 
Each logical qubit is encoded in the decoherence-free sin- 
glet subspace of four spins from the 2x2 plaquette, so 
that it is insensitive to uniform magnetic field fiuctua- 
tions along an arbitrary direction. Besides arbitrary ro- 
tations of single logical qubit achieved by superexchange 
interaction, we provide three different approaches to cou- 
ple the logical qubits from neighboring plaquettes, with 
their properties summarized in Table|T] These approaches 
may also be applied to other quantum systems, such as 
quantum dots or Josephson junction arrays. 
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APPENDIX A: GEOMETRIC PHASE 
APPROACH 

1. Geometric phase approach with bosonic 
particles 

We now justify the claim that the geometric phase tt is 
obtained for 5*2,3 ^Ti/ ^4/ 1 while a trivial geometric phase 
(or 2tt) is obtained for the other three cases (5*2,3 (E" 5*1', 4', 
72,3 (8)51', 4', and T2,3(8)Ti',4/). This evolution implements 
the controUed-phase gate up to a bit-fiip of the logical 
qubit from the left plaquette. 

We start by generalizing the on-site interaction Hamil- 
tonian for site i that governs both the ground and excited 
vibrational levels (denoted as a and b respectively) 



H, 



1. 



1_ 



+ Lu^nt + -Urnt « - 1) + -C/f"^ 



-1) 
(Al) 



U.. 



ab 



where /i^ is the energy off-set, tOi is the vibrational fre- 
quency, U^^ is the on-site interaction strength between 
levels a and (3 for site i. The particle number opera- 
tors are nf = X^o- '^i = J2cr ^i,abi,a, and Ui = 



ab 



+ n^. Given large vibrational frequency uii S> t/, 
we may safely neglect those energy non-conserving terms 



13 




FIG. 12: (Color online) The biased potential for the sites L 
and R. The vibrational ground levels are \L,a) and \R,a), 
with energy difference A. The vibrational excited level for 
the right site is |_R, b) with excitation energy uj. 



{nL,n%) 


3R = {n-k - 1) /2 


jR = (n^ + 1) /2 


(1,0) 
(1,1) 
(1,2) 
(2,1) 
(2,2) 







-W^R 
-2Ur + 2Ur 
-~4Ur + 2Ur 



TABLE II: Energy difference associated with bosonic par- 
ticle tunneling (SEi) for various initial number configura- 
tions [nL,n%) and the final total spin at the right site jr. 
The bias is set to be A = u. The resonance condition 
SEi = is fulfilled for both cases: {nL,n%,jR) = (1,0, 1/2) 
and {nL,n1i,jR) = (1, 1,0) . 



For the biased potential between the two horizontal 



sites (L and R) (as shown in Fig. 12 1, we consider one 
vibrational level for the left site and two levels for the 
right site: 



HL^AnL + UruL {UL - 1) 



(A2) 



Hr = cjnn'n + ^(7^"^ (n?^ - 1) + ^C/j^^n^, (n?, - l) 



1 



(A3) 



where A — — /i^ is the bias in the potential (i.e., 
energy difference between the ground levels for the two 
sites), and Jr is the total spin for the right site (see 
Appendix A 2 for detailed derivation). 

Given quantum numbers (n/,, n^, njj, ifl) , the on-site 
energies are 



El [ul] ^AnL + UriiL [ul - 1) 



(A4) 



Er [nl, 



1, 



where 



(A5) 



jR U 



for 



H > 



jR 



> 



Note that 



1) 
(A6) 

= im- 



plies / [n^, 0,j_R, = 71^/2] = 0. For njj = 1, we have 



n%l:3R 

- i.m 



n„ + l 



In Table |TT] we list energy difference associated with 
bosonic particle tunneling with A — lj. There are two 
possibilities to fulfill the condition of resonant tunnel- 
ing. The first case is {nL,n\,jR) = (1,0,1/2). This 
corresponds to the resonant tunneling between the sites 
(2, 1') in the highlighted upper right panel S'2,3 ® T'i/_4' 
of Fig. [6]^c), yielding a geometric phase tt. The second 
case is lnL,n'^,jR) = (1,1,0). This corresponds to res- 
onant tunneling for both pairs of sites (2, 1') and (3,4') 
as shown in the left panel of Fig. ^c), yielding a triv- 
ial geometric phase 27r. The remaining cases are off- 
resonant \6Ei I ^ t, which yields a trivial geometric phase 
0. Therefore, we obtain a non-trivial geometric phase tt 
only for 5*2 3 ® Ti' 4/ . 



On-site Interaction for Bosonic Particles 



We now derive Eq. (A3 1 from Eq. (Al). The key step 

is to simplify the exchange term a' '^1 a^l a'^i-.o-'^i-.'^' 
using the particle number operators and the total spin 
operator. For clarify, we drop the sub-index i (or R). 

Use the Schwinger representation, we define the spin 
operator J" for the ground vibrational level 



2n^ and/[n?j,l,jfl,= ^ 

Thus the energy difference for the bosonic 
particle tunneling from the left site to the right site is 

6 El [nL,n^, jfl] 

= El [ul] + Er K, 0, ny2] - El {ul -\\-Er K, 1, j^] 
= {A-uj) + Ur {riL - 1) - Vt^ K, 1, Jfl] . (A7) 







Jy = 

















Y' 



(A8) 
(A9) 
(AlO) 
(All) 



with i^Jj = J° {j"" + !)• Similar definition for can 
be introduced for the excited vibrational level. Thus the 
total spin is J = J° -I- J^. 
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We may rewrite the exchange interaction as 



(A12) 



"'U^^:'^4.aK<r + A<y'^^-'^AsK<y (A13) 



4jV + j^ -(/+/) (r+/ + i) 

7i"n^ + - 



(A14) 
(A15) 

(A16) 



Plugging the above expression into Eq. (Al) gives us 
Eq. (|A3l. 



3. Geometric phase approach with fermionic 
particles 

The procedure for the fermionic particles is almost the 
same as that for the bosonic particles, except for the fol- 
lowing three differences. First, the bias of the energy 
off-set needs to be A = + Vj:^ for fermionic particles 
(whereas A = a; for bosonic particles). Second, the geo- 
metric phase TT is obtained from the resonant tunneling 
associated the subspace T2_3(8)5i'4' for fermionic particles 
(whereas it is associated with 5*2^3(8) 11/ 4' for bosonic par- 
ticles). Third, the geometric phase is for the remaining 
cases for fermionic particles (whereas it might be either 
or 27r for bosonic particles). 

For fermionic particles, the on-site interaction Hamil- 
tonian for site i that governs both the ground and excited 
vibrational levels (a and b) is 



where we have safely neglected the energy non-conserving 

terms Y.a,a' ^\,a^\,a'°'h<y<^i,<y' foi' > Uf. 

For the biased potential between the two horizontal 
sites (L and R) as shown in Fig. 12 we have 



ifi = A UL + C/r^^LjnL.i , (A18) 
Hr = un'k + U-R^nX^n%^ + U'j.'n^^n],^^ (A19) 



Given quantum numbers (n^,, n^, n^, ja) , we obtain 
the on-site energy 



(A20) 

En [n?,, n'k, Jr] = cju'r + Ul^S^..^^^ + U'r%,.. (A21) 



El K] = a ni + C/r5„,,2, 
^n-RjR] = 

+ 2- 



(nL,n'k) 


5Ex 


(1,0) 
(1,1) 
(1,2) 
(2,1) 
(2,2) 




1 r rah i rjaa 



TABLE III: Energy difference associated with fermionic par- 
ticle tunneling {5E\) for different initial number configura- 
tions. The bias is set to be A = tj-l-C/Jl''. For (riL, n^) = (1, 1), 
the energy difference is ^If^ for singlet and triplet states, re- 
spectively. The resonance condition SEi = is fulfilled only 
if (ni,n?;j = (l,2). 



where 



Vn 



(3-4jfl)5„»^,i(5„5^^ 



for the spin dependent interaction. If 



(A22) 
I = 1, 



77 = 3 for spin singlet and ry = — 1 for spin triplet states; 
otherwise, r/ = 0. 

The energy difference associated with fermionic parti- 
cle tunneling from the left to the right site is 



5 El [nL,n^, jfl] 
=El [ul] + Er 



-EL[nL-l]-ER[n%ljR] 
(^A-^-ic/]|^(n?j + ,7«?„M«)) 



t/r5n.,2, (A23) 



UJ - 



' Ur, we fulfill the resonance condi- 



By choosing A 

tion SEi = for {nL,n%) = (1,2). In Table [m| we list 



the energy difference associated with fermionic particle 
tunneling for various particle number configurations. 

Here are some remarks on the geometric phase ap- 
proach. It is crucial to have large and sufficiently dif- 
ferent on-site interactions compared with the tunneling 
rate t, because the virtual tunneling process may induce 
higher order systematic errors ~ (t/U) in the accumu- 
lated phase. One may use Feshbach resonances to en- 
hance the on-site interaction and suppress such errors. In 
addition, the tunneling rate t and the energy difference 
A should be as homogeneous as possible, since inhomo- 
geneities St and SA may induce leakage errors out of the 
logical subspace with probability (y) and (^) , re- 
spectively. Optical lattices in flat-bottom traps |3D] may 
efficiently eliminate such inhomogeneity errors. 



APPENDIX B: PERTURBATIVE APPROACH 

1. Effective Hamiltonian with d = J 

For J = d (which might be achieved by placing an ad- 
ditional atom at the center of the plaquette to block its 
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occupancy due to interatomic repulsion |3T]), the singlets 
are degenerated and isolated by an energy gap 4 J from 
the three-fold degenerate S = 1 states. Since they are de- 
generate we can choose any linear combinations of them 
as a basis. For convenience we choose the states |0) and 
|1) as our basis. These two states can be regarded as the 
two components of an effective pseudo-spin 1/2 system. 

In the absence of any coupling the singlet subspace 
of the chain is spanned by product states of the effec- 
tive pseudo-spin states |0),|1) at the plaquettes and all 
four possible configurations are degenerate. A finite J' 
breaks the degeneracy. Since the He Hamiltonian does 
not directly couple the effective pseudo spins, they get 
only coupled through second order virtual processes to 
intermediate high energy states with S > 0. Denoting 
the left and right plaquettes as i and « -I- 1, the effective 
Hamiltonian is given by: 



S,S',q,q' 



2E{0) - E{S) - E{S') 



(Bl) 



Here q{S) labels the number of states within a plaquette 
with total spin S. After some algebra one obtains that 



where a are effective Pauli matrices acting on the !□), | x) 
states, and 



1 


f9J 




8J 


24J 


48 




d 


~ 3 J 


+ 2-1 r + 

d+ J 




f9J 




8J 


^ 2J-d) 


48 




d 


- 3 J 





J 



(B4) 
(B5) 



In Fig. |8]^a) the parameters Xz and 7^ are plotted as a 
function of d/ J. 

The above Hamiltonian is more complex than the one 
derived for the previous J — d case, however if the inter- 
plaquette coupling is smaller than the energy splitting 
between the two singlet states, J'^/J <C AE, it is en- 
ergetically costly to flip an encoded spin and only the 
terms that preserve the total effective magnetization are 
relevant. Consequently, in this regime one can use an 
effective rotating wave approximation which consists on 
neglecting the non-energy-preserving terms (last line in 
Eq. B3l. It leads to a simpler effective Hamiltonian: 



,AE J'^jz 
2 



J 



(B6) 



J 
~7 



'2 



-Si ■ di+i + [Xz - -)o-f(7|+i 



h: 



eff 



3 J 



(B2) APPENDIX C: OPTIMAL QUANTUM CONTROL 
WITH SMOOTH PULSES 



where ct^ are effective Pauli matrices acting on the 
pseudo-spin states |0), |1). H'^^^ is an effective Ising 
hamiltonian and up to single logical qubit rotations it 
can be used to generate cluster states encoded with the 
singlet basis. Since H^di commutes with i?f^{j_|_2, the 
cluster state generation can be performed simultaneously 
in the 2D array of plaquettes. 



Here we present the algorithm used to optimize the 
smooth functions ak- This algorithm is significantly dif- 
ferent from what is the standard approach in optimal 
quantum control based on a Lagrangian formulation [28] , 
and is an extension of the one used by some of the authors 
in j29j. 

Let us consider that we have the Hamiltonian 



2. Effective Hamiltonian with d < J 

For d < J, we shall use the eigen-basis !□) and | x). We 
obtain the effective Hamiltonian using the second order 
perturbation theory: 



K 



K 



H{t;x) = y^Q;fc(t, a;fci. 



fe=i 



k=l 



1=1 



(CI) 

where in our case K — 5 and Ji{t) = sin(Zt7r/T). We 
have defined x = {xki} as the set oi K x L parameters 
that we wish to optimize. The optimization is made such 
that the unitary evolution operator U(T; x) 



/2 r 

J 
J 



(B3) 



^ { I XX \ I ^ \ ^ X-X 



i — U{t;x) = H{t;x)U{t;x) 



(C2) 



with initial condition U{0; x) =1 gets as close as possible 
to a desired unitary gate Ug acting on a subspace {|V'n)} 
of N states of dimension d > N. This can be quantified 
with the fidelity 

1 1 ^ 

/ = -Tr[UlU{T;x)] = - ^(V„|;7tc/(T; x)|V„). 
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To avoid complex numbers one can use either F = Re{/} 
or F = l/p. For simplicity we will derive the algorithm 
using F — Re{/} (but we have presented the results 
with F = l/p as this is not sensible to irrelevant global 
phases). To maximize F we need to compute the deriva- 
tive of F with respect to the parameters x, i.e. dF/dxki- 
If this can be done efficiently, then one can use any of 
the multiple optimization algorithms which compute the 
optimal control. In particular, the derivative and the 
formulas derived below are fed to Matlab's nonlinear op- 
timization toolbox |32j . 

The partial derivatives of F can be expressed 



OF 
dx 



kl 



1 ^ 9 



n=l 



which relates the gradient of to a derivative of the uni- 
tary operator U (T; x) . Using second order perturbation 
theory (see the appendix in [53] for details) the derivative 
of U (T; x) can be expressed as 

^U{t;x) = -iU{t;x) [' dTU{T;x)^^^^^^U{T;x). 
dxki Jo ox 

(C3) 

Hence, we get the following closed formula for the gradi- 
ent of the fidelity 



dxki N ^ 



dr(^„|C/t[/(r)(7(r)t^^C/(r)|^„ 



dx 



kl 



(C4) 

(from hereafter we omit the a;-dependence of U and H 
in order to ease the notation). Though we have a closed 
formula, we still need to perform the integral. To do so, 
we devise an efficient procedure which is based on solv- 
ing three sets of ordinary differential equations (ODEs) . 
First note, that the integral in Eq. (C4| can be trans- 
formed into N X K X L ODEs 



dt 



fnkl{t) 



llm(V;„|C/t[/(T)C/(t)t^C/(i)|V^„ 



(C5) 



with initial conditions fnki(fi) = 0. Thus, we have that 

N 



OF 
dxki 



^Y^UkiiT). 



(C6) 



Then, the algorithm to obtain the gradient of the fi- 
delity, which is fed to Matlab's nonlinear optimization 
toolbox [32j, is given by: 



1. Solve the N ODEs 



^J^\Ut))^Hmn{t)), 



(C7) 



with initial condition |^„(r)) = J7g|'0n) and moving 
backwards in time from T to t. 

2. Solve the 27V ODEs 

AMt)) H{t)\Mt)) (C8a) 
at 

ij^iut)) --^ mmnit)) (C8b) 

with the initial conditions |'0n(O)) = \ipn), and 
|C„(0)) = U{TYUg\tPn) computed before. 

3. Solve the K X L ODEs 

N 



^Im(C„(t)|Ofc|Vn(t)) (C9) 



n=l 



with initial condition ffc/(0) = 0. We have defined 
ffe/ = Sra ^nki and we have used t he ex pression of 
the particular Hamiltonian of Eq. (CI). This step 



can be done simultaneously with step 2 so that 
|'0n(^)) and ICn(i)) need not to be stored. 



4. The derivatives of the fidelity are then given by 
dF 



dx 



kl 



hi{T). 
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